Journal  of  Aircraft 

Vol.  39,  No.  2,  March- April  2002 


Engineering  Notes _ 

ENGINEERING  NOTES  are  short  manuscripts  describing  new  developments  or  important  results  of  a  preliminary  nature.  These  Notes  cannot  exceed  6 
manuscript  pages  and  3  figures;  a  page  of  text  may  be  substituted  for  a  figure  and  vice  versa.  After  informal  review  by  the  editors,  they  may  be  published  within 
a  few  months  of  the  date  of  receipt.  Style  requirements  are  the  same  as  for  regular  contributions  (see  inside  back  cover). 


Effect  of  Strake  Geometry 
and  Centerbody  on  the 
Lift  of  Swept  Wings 

Michael  P.  Schultz*  and  Karen  A.  Flack’ 

United  States  Naval  Academy, 

Annapolis ,  Maryland 21402-5042 

Introduction 

ODERN  fighter/attack  aircraft  rely  on  vortex  flow  to  increase 
lift  at  moderate  to  high  angle  of  attack"  a.  Therefore,  a  great 
deal  of  interest  exists  in  understanding  and  predicting  the  vortical 
flow  structure  generatedby  swept  wings.  One  form  of  vortex  control 
uses  small  geometry  modifications  (fillets)  at  the  strake- wing  junc¬ 
tion  to  alter  the  vortex  dynamics.  Fillets  are  also  possible  candidates 
for  roll  control  devices,  being  deployed  symmetrically  for  enhanced 
lift  and  longitudinalcontrolor  asymmetrically  for  lateral-directional 
control,  as  discussed  by  Hebbar  et  al.1  Excellent  reviews  of  these 
vortex  control  concepts  are  given  by  Lamar,2  Rao  and  Campbell,3 
and  Kern.4 

Various  studies  have  looked  at  the  effect  of  fillets  on  the  interac¬ 
tion  of  the  wing  and  strake  vortices.  A  numerical  study  by  Kern4  on 
a  76/40-deg  swept  wing  indicated  that  linear  and  diamond-shaped 
fillets  could  increase  lift  by  up  to  1 4%  at  low  a  and  up  to  8  %  at  high 
a  with  a  small  improvement  in  lift-to-drag  ratio.  An  experimental 
study  on  the  same  wing  geometry  was  conductedin  a  water  tunnelby 
Hebbar  et  al. 1  The  inclusion  of  fillets  produced  a  more  complicated 
vortex  interaction  caused  by  the  formation  of  additional  vortices  in 
the  fillet  region.  Flow  visualizationindicated  a  delay  in  vortex  burst 
for  diamond- shaped  fillets,  with  no  definitive  trend  in  the  vortex 
response  for  linear  and  parabolic-shaped  fillets.  Another  possible 
benefit  of  fillets  at  the  strake-wing  junction  might  be  a  reduction 
in  tail  buffet.  Impingement  of  burst  vortices  on  the  vertical  tails  of 
modern  fighter/attack  aircraft,  referred  to  as  tail  buffet,  has,  in  some 
cases,  caused  structural  fatigue  and  necessitatedreplacementof  the 
tails.  Ghee  et  al.5  investigated  the  magnitude  and  frequency  of  the 
tail  response  to  vortex  forcing  using  a  matrix  of  pressure  transduc¬ 
ers  on  the  vertical  tails  of  a  76/40-deg  flat-plate,  swept  wing.  The 
parabolic  fillet  created  the  most  tail  buffet  because  this  geometry 
results  in  a  strong  vortex  system.  The  discontinuous  shape  of  the 
diamond-shaped  fillet  results  in  a  number  of  vortices  resulting  in  a 
weaker  vortex  system  that  led  to  the  weakest  tail  buffet. 

The  primary  objective  of  the  present  experiment  was  to  investi¬ 
gate  the  lift  characteristics  of  fillets  on  a  three-dimensional  model 
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and  compare  the  results  to  a  flat-plate  model  with  the  same  plan- 
form.  It  was  also  hoped  to  experimentally  validate  and  quantify  the 
increase  in  lift  as  a  result  of  fillets  on  flat-plate  models  with  the 
same  planformas  the  three-dimensionalmodel.  A  secondary  objec¬ 
tive  was  to  evaluate  the  effect  of  a  simple  centerbody  on  the  lift  of 
flat-plate  models  with  and  without  fillets  as  a  means  for  comparison 
to  the  full  three-dimensionalmodel. 

Experimental  Apparatus 

Wind-Tunnel  Facility 

Tests  were  performed  in  a  subsonic  wind  tunnel  at  the  U.S.  Naval 
Academy.  This  is  a  low-turbulence,  open-return,  closed-jet  facility 
with  a  test  section  79  cm  in  height,  112  cm  in  width,  and  3.05  m  in 
length.  The  tunnel  uses  a  185-kW  electric  motor  controlled  with  a 
variable  frequency  drive  and  delivers  a  maximum  frees  tream  veloc¬ 
ity  of  about  100  ms-1.  Flow  management  in  the  settling  chamber 
consists  of  a  high-aspect-ratio  honeycomb  flow  straightener  and  a 
series  of  fine  mesh  screens.  The  nozzle  contraction  ratio  is  9.6:1. 
The  resulting  freestream  turbulence  intensity  in  the  test  section  is 
0.06%.  The  test  section  is  fitted  with  a  six-component,  sting-type 
force  balance  system.  The  a  of  the  model  can  be  varied  from  —35 
to  +35  deg  using  computer-controlled  stepping  motors.  The  data 
from  the  sting  balance  are  acquired  using  a  Hewlett-Packard  Model 
3852,  13-bit,  data-acquisition  system  connected  to  a  PC. 

Swept-Wing  Models 

A  baseline  swept-wing  model  and  several  derivatives  were  tested 
in  this  investigation.  The  baseline  model  was  a  1 :28  scale,  three- 
dimensional  model  of  a  modern  fighter/attack  aircraft.  The  wing 
span  of  the  model  S  was  45.8  cm  with  a  mean  aerodynamic  chord 
(MAC)  of  14.26  cm  and  a  wing  reference  area  of  592.5  cm2.  It 
was  constructed  of  fiberglass-reinforcedplastic.  Derivatives  of  this 
model  incorporated  linear,  parabolic,  and  diamond-shaped  fillets  at 
the  strake-wing  junctions.  The  fillets  increased  the  wing  area  by 
1%  and  were  similar  to  those  used  by  Kern,4  Hebbar  et  al.,1  and 
Ghee  et  al.5  on  a  76/40-deg  flat-plate,  swept-wing  model.  The  two 
remaining  models  were  flat-plate,  swept-wing  models  (Fig.  la).  The 
first  had  the  same  planform  as  the  baseline  model.  The  second  had 
the  same  planform  as  the  baseline  model  with  the  linear  fillets. 
Linear-shaped  fillets  were  selected  because  they  represent  a  much 
simpler  geometry  to  fabricate  than  the  diamond  or  parabolic  fillets, 
and  the  numericalresults  of  Kern4  showed  that  they  can  significantly 
enhance  lift  at  high  a.  Both  flat-plate  models  were  tested  with  and 
without  the  addition  of  an  ogival  cylinder  centerbody  on  the  top 
surface.  The  centerbody  was  placed  at  a  distance  x /L  =  0  and  0.2 
(Fig.  lb).  These  results  are  referred  to  as  the  forward  centerbody 
and  aft  centerbody  results,  respectively.  All  of  the  flat-plate  models 
used  in  the  experiments  were  7.2  mm  thick  and  were  constructed 
of  biaxial  carbon  fiber  skins  over  a  plywood  core.  The  leading  and 
trailing  edges  of  the  strake,  wing,  and  tail  were  sharp,  20-deg  beveled 
with  a  flat  top  surface. 

Results  and  Discussion 

The  tests  were  conducted  at  a  constant  wind-tunnel  velocity  of 
63.5  ms-1  (Re  =  5.8  x  105  based  on  MAC).  Each  model  was  stati¬ 
cally  pitched  from  — 10  to  +25  deg  in  1-deg  increments  on  a  sting- 
type  six-componentforce-balance  system.  Force-balance  data  were 
collected  at  a  60-Hz  sampling  rate  over  a  1-min  interval.  The  ac¬ 
quired  raw  force  data  were  reduced  in  accordance  with  the  standard 
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Both  Flat— Plate  Models  are  7.2  mm  in  Thickness 

a)  Fillet  details 


No  Centerbody  Aft  Centerbody  Forward  Centerbody 

b )  Centerbody  arrangement 

Fig.  1  Plan  view  of  flat-plate  models. 


procedures  outlined  in  Barlow  et  al.6  This  included  removal  of 
weight  and  deflection  tares,  solid  and  wake  blockage  corrections, 
conversion  to  force  coefficients,  wall  interference  upwash  correc¬ 
tion,  and  conversion  of  force  coefficients  from  body  to  wind  axes. 

Precision  uncertainty  estimates  for  the  resulting  mean  lift  co¬ 
efficients  were  made  using  repeatability  tests.  Ten  replicate  pitch 
sweeps  were  made  on  all  of  the  model  configurations.  The  standard 
error  for  the  mean  lift  coefficient  CL  at  each  a  was  calculated.  To  es¬ 
timate  the  95%  precision  confidence  limits  in  the  mean,  its  standard 
etxor  was  multiplied  by  the  two-tailed  t  value  (t  =  2.262)  for  9  deg 
of  freedom  and  significance  level  of  0.05,  as  given  by  Coleman  and 
Steele.7  These  were  combined  with  the  estimates  of  the  bias  errors 
to  predict  the  overall  95%  confidence  limits  for  CL.  These  were 
estimated  to  be  less  than  ±1.5%  for  all  of  the  models  tested. 

The  lift  coefficient  results  for  several  of  the  configurations  tested 
are  shown  in  Fig.  2.  The  ordinate  in  the  figure  was  chosen  to  be 
Cl~CLo  so  that  the  results  from  the  three-dimensional  and  flat- 
plate  models  could  be  more  easily  compared.  The  lift  coefficient 
at  a  =  0  deg,  CLo,  was  nearly  zero  for  all  of  the  three-dimensional 
cases  tested  and  was  negative  for  all  of  the  flat-plate  cases  because 
of  the  beveled  edge  on  the  suction  side  (see  Fig.  la).  The  center- 
body  and  fillets  did  not  have  a  significant  effect  on  CLo.  The  results 
in  Fig.  2  indicate  that  within  the  experimental  uncertainty  the  lin¬ 
ear  fillet  had  no  significant  effect  on  CL  for  the  three-dimensional 
model  over  the  range  of  a  tested.  Though  not  presented  here,  this 
was  also  the  case  for  the  diamond  and  parabolic  fillets.  These  re¬ 
sults  conflict  with  the  flat-plate  studies  of  Kern4  and  Hebbar  et  al.,1 
which  both  indicated  an  increase  in  lift  with  the  inclusion  of  fil¬ 
lets.  The  three-dimensional  model  results  seem  to  indicate  that 
there  are  distinct  differences  in  the  vortex  interaction  and  break¬ 
down  between  the  flat-plate  and  three-dimensional  model.  The  re¬ 
sults  for  the  flat-plate  model  show  a  significant  increase  in  lift  in 
the  moderate-to-higha  range  with  linear  fillets  as  compared  to  the 
no-fillets  case.  The  increase  in  CL  ranged  from  5.6  to  12.9%  and 
averaged  9.4%  for  a  =  ±10  to  ±20  deg.  The  increase  ranged  from 
2.5  to  6.1%  and  averaged  3.7%  for  a  =  ±20  to  ±25  deg.  These  re¬ 
sults,  although  on  a  different  planform,  show  similar  trends  as  seen 
in  the  numerical  study  of  Kern4  on  a  flat-plate,  76/40-deg  swept 
wing. 

It  was  hypothesized  that  the  differences  in  lift  obtained  for  the 
three-dimensional  model  as  compared  to  the  flat-plate  model  was 
caused  by  a  fuselage  effect.  Results  for  the  flat-plate  model  with 
and  without  a  centerbody  are  also  shown  in  Fig.  2.  The  results  for 
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Fig.  2  Lift  coefficient  results  for  all  model  types  with  and  without 
linear  fillets. 

the  flat-plate  model  with  a  centerbody  show  a  significant  increase  in 
lift  in  the  moderate-to-higha  range  with  linear  fillets  as  compared 
to  the  no- fillets  case.  The  increase  in  CL  ranged  from  7.8  to  12.5% 
and  averaged  9.9%  at  a  =  ±10  to  ±20  deg  for  the  aft  centerbody 
case.  The  increase  ranged  from  2.1  to  7.3%  and  averaged  3.3%  at 
a  =  ±20  to  ±25  deg.  For  the  forward  centerbody  case  the  increase 
in  CL  ranged  from  4.2  to  8.9%  and  averaged  6.0%  at  a  —  ±10  to 
±20  deg.  The  increase  ranged  from  2.2  to  4.2%  and  averaged  3.4% 
at  a  =  ±20  to  ±25  deg.  The  aft  centerbody  had  no  significant  effect 
on  the  lift  for  the  flat-plate  model,  with  or  without  linear  fillets,  as 
these  results  agree  within  their  uncertainty  with  those  taken  without 
a  centerbody.  The  forward  centerbody  results,  with  or  without  linear 
fillets,  indicate  an  increase  in  CL  over  both  the  aft  centerbody  and 
no  centerbody  cases. 

The  results  for  the  three-dimensional  model  and  the  flat-plate 
models  given  in  Fig.  2  show  significant  differences.  Although  this 
is  not  surprising  in  itself,  the  trends  shown  for  the  fillets  vs  no¬ 
fillets  cases  are  not  in  agreement.  As  stated  before,  the  fillets 
increased  the  lift  on  the  flat-plate  models  and  did  not  on  the  three- 
dimensional  model.  This  illustrates  that  caution  should  be  exercised 
when  attempting  to  apply  findings  from  flat-plate  models  to  three- 
dimensional  models  and  full-scale  aircraft.  In  the  present  study  the 
addition  of  a  simple  centerbody  did  not  yield  any  better  agreement 
between  the  model  type.  Further  research  into  size,  geometry,  and 
placement  of  the  centerbody  should  be  undertaken  to  better  address 
the  ability  of  a  flat-plate  model  to  predict  the  trends  in  performance 
of  a  three-dimensional  model. 

Conclusions 

The  following  conclusions  can  be  made  from  the  present  experi¬ 
mental  investigation: 

1)  Linear  fillets  placed  at  the  strake-wing  junction  on  a  flat-plate, 
swept-wing  model  led  to  a  significant  increase  in  CL  at  moderate- 
to-higha  (±10  to  ±25  deg). 

2)  Linear,  parabolic,  and  diamond  fillets  at  the  strake-wing  junc¬ 
tion  on  the  three-dimensionalmodel  had  little  effect  on  CL  over  the 
entire  a  range  (— 10  to  ±25  deg). 
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3)  A  centerbody  on  the  flat-plate  models  had  little  effect  on  the 
lift  coefficient  and  did  not  produce  better  agreement  with  the  three- 
dimensional  model  than  the  flat-plate  model  alone. 

4)  Caution  should  be  exercised  when  attempting  to  apply  re¬ 
sults  and  overall  trends  from  flat-plate,  swept- wing  models  to  three- 
dimensional  models  and  full-scale  aircraft. 
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Introduction 

HE  need  for  addressing  optimization  problems  that  are  charac¬ 
terized  by  the  presence  of  a  large  number  of  design  variables, 
complex  constraints,  and  discrete  design  parameter  values  exists  in 
many  fields  including  engineering  design.  A  variety  of  local  and 
global  optimization  algorithms  have  been  developed  for  addressing 
such  problems.  Besides  deterministic  methods,  stochastic  methods 
such  as  genetic  algorithm  (GA)  and  simulated  annealing  (SA)  al¬ 
gorithm  have  recently  found  applications  in  many  practical  engi¬ 
neering  design  optimization  problems.  These  algorithms  are  easily 
implemented  in  robust  computer  codes  as  compared  with  deter¬ 
ministic  methods  because  they  do  not  depend  on  direct  gradient 
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information,  which  most  deterministic  methods  do.  However,  SA 
and  GA  methods  require  a  large  number  of  function  evaluationsand 
relatively  longer  computation  time  than  deterministic  methods,  es¬ 
pecially  in  the  case  of  complex  design  problems.  Although  the  use 
of  parallel  GA  and  parallel  S  A  as  outlinedin  Wang  and  Damodaran1 
offers  a  way  to  reduce  the  large  computational  time,  an  attractive 
alternative  to  SA  and  GA  could  be  the  simultaneous  perturbation 
stochastic  approximation  (SPSA)  method  described  in  Spall.2  The 
SPSA  method  has  been  applied  to  numerous  difficult  multivariate 
optimization  problems  in  many  diverse  areas  such  as  statistical  pa¬ 
rameter  estimation,  feedback  control,  simulation-based  optimiza¬ 
tion,  signal  and  image  processing,  and  experimental  design.  The 
essential  feature  of  SPSA,  which  accounts  for  its  power  and  relative 
ease  of  implementation,  is  the  underlying  gradient  approximation, 
which  requires  only  two  measurements  of  the  objective  function  re¬ 
gardless  of  the  dimensions  of  the  optimizationproblem.This  feature 
allows  for  a  significant  decrease  in  the  cost  of  optimization,  espe¬ 
cially  for  problems  with  a  large  number  of  variables  to  be  optimized. 

The  aim  of  this  Note  is  to  compare  performance  of  SPSA  with 
SA  and  GA  and  to  explore  any  advantages  that  SPSA  might  offer 
to  overcome  the  large  computational  efforts  of  SA  and  GA  when 
applied  to  wing-designproblems.  These  methods  are  briefly  outlined 
following  the  statement  of  the  wing-design  optimization  problem, 
which  will  form  the  application  problem  to  assess  and  compare  the 
performance  of  SPSA  in  relation  to  SA  and  GA. 

Wing-Design  Problem 

The  application  concerns  the  design  of  wing  shape  such  that  the 
aerodynamic  efficiency  of  the  wing  or  the  lift  L  to  drag  D  ratio 
reaches  a  maximum  value  during  cruise  with  the  wing  weight  act¬ 
ing  as  a  constraint,  that  is,  the  goal  is  to  determine  the  wing  geom¬ 
etry  by  either  minimizing  D/L  or  maximizing  L/D  with  the  wing 
weight  as  a  constraint.  The  D/L  ratio  can  be  formulated  in  detail 
using  the  analytic  formulas  for  aerodynamic  analysis  as  defined  in 
Raymer.3  The  lift  L  is  defined  as  L  =  CLqS,  where  q  =  \pV2 is 
the  dynamic  pressure,  p  is  the  density  of  air,  V  is  the  flight  speed, 
CL  =  CLa a  is  the  the  lift  coefficient  where  a  is  the  angle  of  attack 
andCLa  =  2jtAr/{2  +  s/[A  +  (AR/)/q)2(l  +  tan2  X  /yS2)] }  is  the  lift 
curve  slope.  In  the  expression  for  lift  curve  slope,  Ar (=  b2 /S)  is 
the  wing  aspect  ratio,  where  b  is  the  wing  span,  X  is  the  wing 
sweep  angle,  q  (value  of  which  lies  in  the  range  0.95-1.0)  is  the 
airfoil  efficiency  factor,  /?  =  1  —  M2  is  the  compressibility  fac¬ 
tor,  and  M  is  the  Mach  number.  The  total  drag  is  defined  as 
D  =  CDqS,  where  the  total  drag  coefficient  is  CD  =  CDl  +  CD0, 
which  consists  of  the  induced  drag  coefficient  CDi  =  C2L/(n ARe) 
and  the  zero-lift  drag  coefficient  CD0  —  CfFQ.  In  these  expressions 
e  =4.61(1  —  0.045A°-68)(cosk)015  —  3.1  is  the  wing  planform 
efficiency  factor,  Cf  =  0.455/[(log10  Re)258(l  +  0.144M2)0-65] 
is  the  surface  skin-friction  coefficient,  which  is  a  func¬ 
tion  of  the  Reynolds  number  Re,  F  =  {1  +  [6.6/(x /c)m](t  /c)  + 
100(//c)4}[1.34M018(cosA)0-28]  in  which  t/c  is  the  airfoil 
thickness-to-chord  ratio,  ( x/c)m  is  the  chord- wise  location  of  the 
maximum  thickness-to-chord  ratio,  taken  as  0.3  in  the  present 
study,  and  Q  is  a  factor  accounting  for  interference  effects  on 
drag  taken  as  1 .0  in  the  present  study.  The  weight  of  the  wing  (in 
pounds)  is  l¥wing  =0.0106(lVdgA,)a5  50-622A«75(//c)-0-4(cosk)-1, 
where  VP^is  the  design  gross  weight  in  pounds  and  Nz  is  the  ultimate 
load  factor,  which  is  assumed  to  be  13.5  for  subsonic  flow. 

The  design  variables  for  the  wing  design  optimization,  that  is, 
a,  b,  c,  X,  and  IPwing,  represent  the  angle  of  attack,  wing  span,  mean 
aerodynamicchord,  sweep  angle,  and  wing  weight,  respectively. The 
objective  function  to  be  optimized  is  F(X)  —  D/L  and  is  defined 
as  follows: 

Minimize  F(X)  (1) 

subjectto  six  constraintson  the  design  variables  defined  as  follows: 

1.0  deg  <aS  10.0  deg,  10.0  <  b  <  50.0 

3.5  £  c  <  10.0,  O.Odeg  <  X  <  35.0deg 

0.5  <  Ar  <  15.0,  Wwing  <  2473  (lb)  (2) 


